Cotunneling Spectroscopy in Few-Electron Quantum Dots 
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Few-electron quantum dots are investigated in the regime of strong tunneling to the leads. Inelastic 
cotunneling is used to measure the two-electron singlet-triplet splitting above and below a magnetic 
field driven singlet-triplet transition. Evidence for a non-equilibrium two-electron singlet-triplet 
Kondo effect is presented. Cotunneling allows orbital correlations and parameters characterizing 
entanglement of the two-electron singlet ground state to be extracted from dc transport. 



Transport studies of few-electron quantum dots have 
proven to be a rich laboratory for investigating the ener- 
getics of electrons in artificial atoms 0, S Q as well as 
related spin effects, including ground- state spin transi- 
tio ns P^IM IEB. spin lifetimes 8] and Kondo effects 
0, lid The interplay of electron-electron interac- 

tions, electron spin, and coupling to a Fermi sea makes 
transport in the few-electron regime a subtle problem in 
many-body physics 0, E Q E H E| • Of particular 
importance is the two-electron case ("quantum dot he- 
lium") ^3 since this is a paradigm for the preparation 
of entangled electronic states 18j, and in double quan- 
tum dots is the basis of a quantum gate proposal JJJ. 

In this Letter, we present a detailed experimental in- 
vestigation of cotunneling through quantum dots con- 
taining one, two, and three electrons. Measurements 
of inelastic cotunneling are used to extract the singlet- 
triplet (ST) splitting across the two-electron ST tran- 
sition. Evidence of a non-equilibrium ST Kondo effect 
for two electrons is presented. Cotunneling and Kondo 
effects are used to determine the g-factor for magnetic 
fields along different directions in the plane of the 2D elec- 
tron gas (2DEG), giving isotropic g-factors close to the 
bulk GaAs value. Using both cotunneling and sequential 
tunneling data, we extract quantum correlations of the 
two-electron singlet ground state, allowing the degree of 
spatially separated entanglement to be measured. 

Previous transport studies of few-electron quantum 
dots have identified the ST ground state transition for 
two electrons S 0? S 01 ^s well as for larger electron 
numbers [23, E| • Inelastic cotunneling was recently in- 
vestigated in few-electron vertical structures in Ref. . 
These authors demonstrated that inelastic cotunneling 
provides a direct and sensitive measure of excited state 
energies. Here, we use this fact to measure the ST split- 
ting, J, across the ST transition (for both negative and 
positive J), and for the first time extract two-electron 
ground state wave function correlations from cotunnel- 
ing. 

Transport through the ST transition has been studied 
theoretically 12], with a prediction of enhanced Kondo 
correlations at the ST crossing 13]. Effects of lifting 
spin degeneracy of the triplet have also been theoreti- 
cally investigated |1^]. For the degenerate triplet case. 



a characteristic asymmetric peak in conductance at the 
ST crossing has been predicted 0, EB- This predicted 
asymmetric peak is observed in the present experiment. 
Previous measurements of ST Kondo effects 0, 0, |2^ 
in dots have not treated the two-electron case. 

Measurements were carried out on two similar lateral 
quantum dots formed by Ti/Au depletion gates on the 
surface of a GaAs/Alo.sGao.rAs heterostructure 105 nm 
above the 2DEG layer (Fig. 5, inset). The two devices 
showed similar results; most data are from one of the 
dots, except those in Fig. 5. The dilution refrigerator 
base electron temperature was T^i = 45 ± 5 mK, mea- 
sured from Coulomb blockade peak widths. Differential 
conductance g = dl/dVsD was measured with typical ac 
excitations of 5 /jY. 
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FIG. 1: (a) Differential conductance g (log color scale) as a func- 
tion of source-drain bias Vsd and gate voltage Vq at B± = 0, 
at base electron temperature T^i = 45 mK. Numbers through 4 
are number of electrons in the dot. White vertical lines identify 
the locations for data shown in (c) and (d). (b) Same as (a), at 
B± = 1 T. (c) Differential conductance through the N = 2 dia- 
mond showing step with overshoot at Vsd = J(B±)/e at Bj^ = 
and 1 T. (d) Differential conductance through the = 3 diamond 
showing Kondo peak at Vsd =0 for B = 0, split by B± = 1 T. 
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FIG. 2: (a) Differential conductance g as a, function of Vsd in 
the = 1 diamond (Vq = 0.1 V) for in-plane fields Bx = 
0,0.4,0.6,0.8,1 T, (top to bottom, curves offset). Dashed grey 
lines are guides to the eye showing the cotunneling gap. (b) g(VsD) 
shows a zero-bias peak in the = 3 valley (Vq = 0.42 V) that splits 
in an in-plane field By = 0,0.25,0.45,0.7,0.95 T (top to bottom, 
curves offset). (c,d) splitting energies (see text) versus magnetic 
field as in (a,b) with linear fits. Insets: angular dependence of 
the g-factor in the plane of the 2DEG indicating isotropic behav- 
ior. Dashed circles show direction-averaged g-factors. Directions 
X and Y in the plane are arbitrary. 
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FIG. 3: (a) Differential conductance g (log color scale) as a func- 
tion of Vs D and Vg for B = T in the vicinity of the = 1 ^ 2 
transition, (b) g(ysD^B±^ at Vq = 0.164 V (vertical white line in 
(a)) 27J shows the perpendicular field dependence of the singlet- 
triplet gap. (c) Cuts showing g(VsD) at the positions of the verti- 
cal lines in (b), marked A, B, C, D. (d) Horizontal cut (E in (b)) 
showing g(B_\_) at zero bias. Note the asymmetric peak in g at the 
singlet-triplet transition. 



Figures l(a,b) provide an overview of transport spec- 
troscopy data. Diamond patterns of high conductance 
correspond to gate voltages Vg where the ground state 
of the dot ahgns with the chemical potential of either the 
source or drain, allowing sequential tunneling through 
the dot 0, . Transport is absent at more negative 
gate voltages, indicating the absolute occupancy of the 
dot {N = to 4). Conductance features that vanish 
below a finite source-drain voltage \Vsd\ = involve 
transport through an excited state at energy A above 
the ground state. An example of the latter is the nearly 
horizontal band running through the center of the N = 2 
diamond. Beyond this band transport through the ex- 
cited triplet channel of the N = 2 dot becomes allowed, 
as discussed below. 

Inside the diamonds, sequential tunneling is Coulomb 
blockaded and transport requires higher order (cotunnel- 
ing) processes 22, 25]. Elastic cotunneling leaves the en- 
ergy of the dot unchanged; inelastic cotunneling, which 
leaves the dot in an excited state, requires energy sup- 
plied by the source-drain bias. The inelastic mechanism 
becomes active above a threshold Vs d and is independent 
of Vb. 

We first discuss the one-electron regime. A conduc- 
tion threshold within the = 1 diamond [Figs. l(a,b)] 
emerges from the crossing of ground-state and excited- 
state sequential tunneling lines These features cor- 
respond to the onset of inelastic cotunneling through the 
first orbital excited state lying Ai ~ 1.2(1.0) meV above 
the ground state for a field B± = 0(1) T perpendicular to 
the 2DEG. Measurements with magnetic fields up to 1 T 



along different directions in the plane of the 2DEG show 
inelastic cotunneling through Zeeman split one-electron 
states [Fig. 2(a)]. Measurement of Zeeman energies via 
cotunneling was established in Ref. [2^. The cotunnel- 
ing gap Acot — extracted by taking half the peak splitting 
of dg/dVsD — is shown in Fig. 2(c) for one of the field di- 
rections. The g-factors are extracted from a linear fit to 
Acot(^) and are found to be isotropic within experimen- 
tal error, giving a value of (g) = 0.40 ±0.03 averaged over 
the measured field directions. This is close to the bulk 
GaAs value and consistent with previous (few-electron) 
experiments 1] . 

For A^ = 3, a zero-bias conductance peak, presumably 
due to the Kondo effect 9], splits in both perpendicular 
[Fig. 1(d)] and in-plane [Fig. 2(b)] magnetic fields. The 
splitting Ak due to in-plane field — taken as half the dis- 
tance between maxima of the split peaks [indicated in 
Fig. 2(b)] — is shown in Fig. 2(d) along with a best fitting 
line. Slopes from the fits do not depend on direction in 
the plane, and give (g) = 0.43 ± 0.03, consistent with 
the one-electron cotunneling data [Fig. 2(c)]. Note that 
unlike the cotunneling data, the Kondo data does not 
extrapolate to Ak(0) = 0, as also reported in previous 
experiments |26|]. The threshold in-plane field Bk for 
the appearance of Kondo peak splitting gives an esti- 
mate of the Kondo temperature {gjiBBx ^ ^bTk) of 
Tk - 150 mK 26]. 

A detailed view of two-electron transport is shown in 
Fig. 3(a). The nearly horizontal band running through 
the N = 2 diamond [see also Fig. l(a,b)] corresponds to 
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FIG. 4: (a) First and second one-electron excited state energies 
Ai and A2, measured from sequential tunneling along with fits to 
a 2D anisotropic harmonic oscillator model with huja = 1-2 meV 
and hujij = 3.3 meV (see text), (b) Singlet-triplet splitting J as a 
function of magnetic field B± . (c) Dependence of J on gate voltage 
Vg for various B± as indicated, (d) Average slopes AJ(Vg)/ AVq 
from (c) as a function of magnetic field Bj^ , showing strong reduc- 
tion of gate voltage dependence of J at large B±. 



the onset of inelastic cotunneling through the triplet ex- 
cited state, which becomes active for \Vsd\ > J/e. The 
inelastic cotunneling edges align with the triplet excited 
state lines seen in sequential tunneling outside the dia- 
mond, as expected 22]. We use this cotunneling feature 
to measure the ST splitting J. The zero-field value mea- 
sured here, J{B = 0) ~ 0.2 meV, is much less than 
the N = 1 orbital level spacing due to strong interac- 
tions, consistent with theoretical estimates [l^ and pre- 
vious measurements 5]. A surprising zero-bias conduc- 
tance peak in the middle of the cotunneling gap, visible 
in Fig. 3(a) in the range 0.12 V< Vg < 0.15 V is not 
understood. 

Perpendicular field dependence of the ST splitting 
J{B±) is investigated by plotting g along a cut through 
the N = 2 valley as a function of B± [Fig. 3(b)]. Near 
5x = 5* - 1.3 T the ST gap closes and then re- 
opens at larger fields. We interpret this as a ST cross- 
ing where the triplet state becomes the ground state for 
B± > B* [Fig. 4(b)]. We note that in-plane fields up to 
1 T cause no observable change in the two electron spec- 
trum. We also find that J depends on t he g ate voltage 
Vg [Fig. 4(c)], as observed previously jsl. l20j . though at 
larger fields this dependence becomes significantly weaker 
[Fig. 4(d)]. The zero-bias conductance within the N = 2 
diamond as a function of field shows a large, asymmetric 
peak at B± = 5* [Fig. 3(d)], consistent with predictions 
for elastic cotunneling at the ST crossing W} (see also 

His]). 

Before turning to wave function correlations, we first 
extract some useful information about the dot shape from 
the = 1 excitation spectrum. Transport spectra for the 
N = ^ 1 transition, extracted from plots like Fig. 1(a) 
in the region between the N = and N = 1 diamonds, 
give first (second) excited state energies lying A 1(2) above 
the ground state. We find A2 ~ 2Ai, indicating roughly 



FIG. 5: Differential conductance as a function of Vsd for tem- 
peratures Tel = 45, 140, 170, 240, 280, 330, 380, 420, 570 mK (top to 
bottom) showing overshoot at Vsd ~ J/^- Inset: peak conduc- 
tance as a function of temperature with best-fit log{T) dependence 
(line). 



harmonic confinement. Dependencies of Ai(2) on Per- 
pendicular field are well described by a 2D anisotropic 
harmonic oscillator model [2^. From zero-field data, we 
extract fiuoa = 1.2 meV where a{h) is along the larger 
(smaller) dimension of the dot; the energy scale for the 
smaller direction is found by fitting the field dependence 
of Ai(5^), which gives Tiuo}) = 3.3 meV 28j. As a check 
of these values, good agreement between experimental 
and predicted values for /S.2{B±_) is found [Fig. 4(a)]. We 
conclude that the dot potential is spatially elongated by 
a factor of ~ 1.6 = ^Juj^fuja- 

We note that for strong coupling of the dot to the leads, 
the onset of inelastic cotunneling at Vsd = J/e shows 
considerable overshoot, as seen in Fig. 5 (measured in a 
device similar to the one discussed above, with larger ST 
splitting, J(0) ~ 0.57 meV). The temperature depen- 
dence of the maximum overshoot is shown in the inset 
of Fig. 5 along with a line indicating a Kondo-inspired 
log(T) dependence jlO, 111123]. The FWHM of the cor- 
responding positive peak in dg/dVsD is proportional to T 
at high temperatures and saturates at T ~ 80 mK, giving 
an estimate of Tk for this device. However, a quantita- 
tive theory of nonequilibrium ST Kondo effect would be 
needed to further analyze these data. 

Finally, we investigate correlations in the two-electron 
wave function following the analysis of Ref. [l^ . We note 
that Ref. specifically considers a two-electron double 
quantum dot; we anticipate that the elongated shape of 
our single dot will lead to a spatially separated charge ar- 
rangement for A/" = 2, not unlike a double dot in the limit 
of strong interdot coupling. Selecting basis states appro- 
priate for a double dot but applicable here as well — i.e., 
symmetric (|+)) and antisymmetric (| — )) states along 
the long axis of the dot — we identify |+) and | — ) with the 
orbital ground and first excited states of the one-electron 
dot. Because of electron-electron interactions, the N = 2 
ground-state singlet generally comprises an admixture of 
the one-electron ground and excited orbital states. The 
amount of admixed excited state | — ) is parameterized 
by (0 < < 1), the so-called interaction parameter. 
Knowing allows two other important quantities to be 
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extracted: the double occupancy, D = (1 — ^)^/2(1 + ^^), 
and the concurrence 29], c = 2(/)/(l + (/)^), which respec- 
tively parameterize correlations and entanglement of the 
two-electron singlet ground state 16]. 

To extract <p from elastic cotunneling data, one also 
needs to know the charging energy for adding the second 
electron, the operating position within the N = 2 dia- 
mond, and the couplings to each lead for both the singlet 

S T 

and the triplet, • fields well below the ST tran- 
sition, these r's can be estimated from excitation spec- 
tra at the = 1 ^ 2 transition by fitting a thermally 
broadened Lorentzian to the tunneling lines |30|. Upon 
inserting these quantities into Eqs. 8 and 10 of Ref. [Tg, 
we find ^ ~ 0.5 ± 0.1, indicating that the N = 2 ground- 
state singlet contains a significant admixture of the ex- 
cited one-electron orbital state due to electron-electron 
interactions. We emphasize that this method does not 
rely explicitly on a double dot interpretation 31]. From 
this value of (j) we extract a concurrence of c ~ 0.8 for 
the two-electron singlet. This is close to the maximum 
concurrence value c = 1, which characterizes a pair of 
singlet-correlated electrons in fully non-overlapping or- 



bital states. 

Two alternative methods for estimating (j) give consis- 
tent results with the cotunneling method. First, one roay 
adapt the formula (j) = y^l + (u/Uf - U/U from [ig 
by associating the measured Ai with the tunnel split- 
ting 2t of the two lowest noninteracting single-particle 
states, and the charging energy to add the second elec- 
tron with U. The second alternative method uses the 
size of the elastic cotunneling step at the ST transi- 
tion [see Fig. 3(d)] which is shown to be related to (j) 
in 16]. It is notable that all three methods allow the 
concurrence, a measure of "useful" (i.e., spatially sepa- 
rated) two-particle entanglement, to be extracted from a 
dc transport measurement. 
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